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Abstract 
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The two-loop electroweak bosonic correction to the muon lifetime and the My/-Mz inter- 
dependence is computed using analytical methods of asymptotic expansion. Combined with 
previous calculations this completes the full two-loop correction to Ar in the Standard Model. 
The shift to the prediction of VF-boson mass due to two-loop bosonic corrections does not 

CN ' exceed 1 MeV for the range of the Higgs boson mass from 100 to 1000 GeV. 
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O The Fermi constant Gp plays an important role in the precision tests of the Standard Model. 

q ■ Theoretically Gp can be related to other precision observables: the electroweak coupling constant a 

' and the masses of electroweak gauge bosons Mz and M\y- Other parameters enter this expression 
P-T 

i ' through quantum corrections. Usually one inverts this relation in order to predict Mw through Mz 
Qh! 

D . which is measured much more accurate. This Mz — My/ interdependence can be then confronted 

, 

!• ■ with experimental value M^ p . The current error (39 MeV) of M^ p will be drastically reduced at 
^ ' future colliders. In fact, at LHC the experimental error can be reduced to 15 MeV [1] and at Linear 
\ Collider even down to 6 MeV [2] . Therefore much efforts have been spent to reduce the error of the 
theoretical prediction. 

Fermi constant is defined as the coupling constant in the low energy four fermion effective 
Lagrangian describing the decay of the muon 

G 

Cf = £qed + -j= [jyy°(i - 7s)/i] [er«(l - 7sK] (1) 
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where e and \i are electron and muon fields, v e and are the corresponding neutrinos and Gp is 
the Fermi constant. From the Lagrangian (1) one gets the following value for the muon lifetime 



2™5 / 9 



where the factor Ag describes all the quantum corrections in the low energy effective theory (i.e. 
QED corrections). At one- loop order these corrections have been computed a long time ago [3]. 
Recently also the two-loop result for Aq has been obtained [4]. Taking it into account, the error of 
Gp is nowadays dominated by the experimental error of r M measurement 5 . 

In order to relate Gp to My/ one can use the matching condition between effective theory (1) 
and the Standard Model, which requires that the value of r M does not depend on whether it is 
evaluated in the Fermi theory or in the full Standard Model. This brings us to the relation 

Gp e 2 
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Ar), (3) 



where e is the electric charge, sw = sin 9 W is the SM mixing angle and M w is the PF-boson 
mass. The factor e 2 /8s 2 v m 2 v corresponds to the matching at the three level, while the quantity 
Ar = Ar 1_loop + Ar 2_loop + ... is defined to absorb quantum corrections coming from the matching 
procedure at higher loop orders. At the one- loop level Ar was first evaluated in [5]. Further 
improvement of Ar was achieved in [6] where the leading and subleading large t-quark corrections 
were computed at 0(a 2 ) order. Other two-loop fermionic corrections were considered in a series of 
papers [7]. 

The aim of this work is to compute the remaining 0(a 2 ) pure bosonic contribution in order to 
complete the full 0(a 2 ) approximation to Ar. As this work was in progress Ref. [8] has appeared 
where this missing contribution has been computed numerically using the integral representation 
for the massive two-point functions [9]. 

Our approach is complementary to that of [8] . We use the methods of asymptotic expansion [10] 
in two different regimes: 1) large Higgs mass expansion and 2) expansion in difference of masses 
Mjj — Mf . In addition we consider also s 2 ^ = sin 2 9 W as a small parameter in the spirit of [11]. 

The details of our calculation will be presented in a separate publication [12] and here we only 
sketch the most important moments. 

• Gf is nothing but the Wilson coefficient function of the corresponding operator in the low 
energy four fermion effective theory. Therefore it depends only on the "hard" physics and is 
insensitive to low energy field modes in loops. The simplest way to compute such a quantity 



5 The present value is Gf = 1.16637(1) x 10 5 GeV 2 and possible future experiments could reduce the error by 
an order of magnitude. 
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is the use of the factorization theorem (if it exists), which allows one to separate "hard" and 
"soft" physics 6 . This theorem for the construction of low energy effective theories is known 
for a long time (see e.g. [13]) and is based on the euclidean structure of large mass expansion. 

• Using the factorization theorem the problem is reduced to the evaluation of two-loop bubble 
diagrams. 

• All the calculations have been performed in a general gauge with three arbitrary gauge 
parameters. By explicit calculation we checked that contribution to Ar is gauge invariant. 

• In order to have explicit gauge invariant expressions also at the level of bare quantities we 
explicitely include tadpole diagrams. We have checked that applying this prescription all 
the bare parameters and all the counterterms for the masses and the coupling constant are 
explicitely gauge invariant. 

• The renormalization was performed in two different schemes: MS and on-shell scheme. For 
the transformation from one scheme to another one needs the results of [11] where the trans- 
formation formulae for the gauge boson masses were given to two loops. 

To obtain the two-loop bosonic contribution to Ar a number of two-loop diagrams of the order of 
10 4 should be evaluated. In our calculation the computer algebra system FORM [14] was used. In 
order to obtain the FORM readable input we make use of the input generator DIANA [15]. 
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Figure 1: Dependence of Ar^ s on Higgs mass M# (see explanation in the text). Parameter choice 
for the pole masses: M z = 91.1876 GeV and M w = 80.423 GeV. 

6 In our case by "soft" we understand energies and momenta of order of muon mass m M , while "hard" means scales 
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Let us now present the results of the calculation. In Fig. 1 the results of expansions are plotted 
as a function of the on-shell Higgs boson mass Mh starting from the experimental low bound 114 
GeV [16] up to 1 TeV. First we expand in the limit M H — > M z , e.g. in the difference (M H — M§). 
By summing this series we have found that with 6-10 coefficients it works well for the values of 
Mh down to ~ 125 GeV. In order to accelerate the convergence of this series we apply Pade 
approximation w.r.t. the variable y = (M H — M z )/M z . The left bold curve in Fig. 1 represents the 
[3/3] Pade approximant which behaves smooth down to ~ 170 GeV. The other bold curve is a sum 
of the first 4 terms of large Higgs mass expansion. This series appeares to be nonalternating and 
the application of Pade did not show any considerable improvement. In the intermediate region 
155 GeV < Mh < 350 GeV both approximations fail to converge. (This is due to the presence of 
singularities in M^-plane at |M#| = 2Mw or \Mh\ = ZMz- Such singularity appears for example 
in diagram of Fig. 2. This manifests itself as a divergence of the series as Mh approaches 2M Z -) 
Therefore, in order to get a result for Higgs mass values in the region 155-350 GeV we made simple 
polynomial interpolation. 
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Figure 2: Example of a diagram that has singularity near Mh = 2M\y- 



Our comparison with the numerical curve from Ref. [8] shows a perfect agreement between the 
two calculations. The difference for 155 GeV < M H < 350 GeV does not exceed 5% while for M H 
below 155 GeV and above 350 GeV it is much smaller. 

Let us now turn to the Mz — My/ interdependence. By inverting Eq. (3) we can find how much 
does the computed Ar££, contribute to the shift of the value of the VF-boson mass. This contribution 
is plotted in Fig. 3. For the broad region of M H from 100 to 1000 GeV this correction does not 
exceed 1 MeV and appears to be insignificant. Therefore, we can see that the uncertainty in M w 
prediction at the moment is dominated by the uncertainty in the experimental value of t-quark pole 
mass. 

In conclusion, the bosonic contributions to Ar has been computed using methods of asymptotic 
expansions and the low energy factorization theorem. The influence of this correction on the theo- 
retical prediction of M w is found. The gauge invariance of Ar^ s has been proved explicitly. By the 
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Figure 3: Shift AM W of the predicted mass of W^-boson due to two-loop bosonic correction to Ar. 

inclusion of Higgs-tadpole diagrams it has been checked that all mass- and charge- counterterms 
are gauge invariant. 
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